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00 . ABSTRACT 



We continue our work on developing techniques for studying turbulence with spectroscopic data. 
We show that Doppler-broadened absorption spectral lines, in particularly, saturated absorption lines, 
^ ^ ' can be used within the framework of the earlier-introduced technique termed the Velocity Coordinate 

^ \ spectrum (VCS). The VCS relates the statistics of fluctuations along the velocity coordinate to the 

. statistics of turbulence, thus it does not require spatial coverage by sampling directions in the plane of 

the sky. We consider lines with different degree of absorption and show that for lines of optical depth 
OO \ less than one, our earlier treatment of the VCS developed for spectral emission lines is applicable, if 

the optical depth is used instead of intensity. This amounts to correlating the logarithms of absorbed 
intensities. For larger optical depths and saturated absorption lines, we show, that the amount of 
information that one can use is, inevitably, limited by noise. In practical terms, this means that 
only wings of the line are available for the analysis. In terms of the VCS formalism, this results in 
O ' introducing an additional window, which size decreases with the increase of the optical depth. As a 

result, strongly saturated absorption lines carry the information only about the small scale turbulence. 
^ . Nevertheless, the contrast of the fluctuations corresponding to the small scale turbulence increases with 

^ ' the increase of the optical depth, which provides advantages for studying turbulence combining lines 

with different optical depths. We show that, eventually, at very large optical depths the Lorentzian 
profile of the line gets important and extracting information on velocity turbulence, gets impossible. 
^ ' Combining different absorption lines one can tomography turbulence in the interstellar gas in all its 

complexity. 

' Subject headings: turbulence - ISM: general, structure - MHD - radio Hnes: ISM. 



1. INTRODUCTION 

Turbulence is a key element of the dynamics of astrophysical fluids, including those of interstellar medium (ISM), 



\ clusters of galaxies and circumstellar regions. The realization of the importance of turbulence induces sweeping changes, 
• • . for instance, in the paradigm of ISM. It became clear, for instance, that turbulence affects substantially star formation, 
. ^ mixing of gas, transfer of heat. Observationally it is known that the ISM is turbulent on scales ranging from AUs to 
Jv>( , kpc (see Armstrong et al 1995, Elmegreen & Scalo 2004), with an embedded magnetic field that influences almost all 
. of its properties. 

' The issue of quantitative descriptors that can characterize turbulence is not a trivial one (see discussion in Lazarian 
1999 and ref. therein). One of the most widely used measures is the turbulence spectrum, which describes the 
distribution of turbulent fluctuations over scales. For instance, the famous Kolmogorov model of incompressible 
turbulence predicts that the difference in velocities at different points in turbulent fluid increases on average with the 
separation between points as a cubic root of the separation, i.e. \5v\ ^ V-^^. In terms of direction-averaged energy 
spectrum this gives the famous Kolmogorov scaling E{k) ~ 47rfc2p(k) - fc^/s^ ^here P(k) is a 3D energy spectrum 
defined as the Fourier transform of the correlation function of velocity fluctuations f(r) = ((5t;(x)(5w(x -I- r)). Note that 
in this paper we use (...) to denote averaging procedure. 

Quantitative measures of turbulence, in particular, turbulence spectrum, became important recently also due to 
advances in the theory of MHD turbulence. As we know, astrophysical fluids are magnetized, which makes one believe 
that the correspondence should exist between astrophysical turbulence and MHD models of the phenomenon (see 
Vazquez-Seniadeni et al. 2000, Mac Low & Klessen 2004, Bellesteros-Paredes et al. 2007, McKee & Ostriker 2007 and 
ref. therein). 

In fact, without observational testing, the application of theory of MHD turbulence to astrophysics could always be 
suspect. Indeed, from the point of view of fluid mechanics astrophysical turbulence is characterized by huge Reynolds 
numbers, i?e, which is the inverse ratio of the eddy turnover time of a parcel of gas to the time required for viscous 
forces to slow it appreciably. For Re 3> 100 we expect gas to be turbulent and this is exactly what we observe in HI 
(for HI Re ^ 10*). In fact, very high astrophysical Re and its magnetic counterpart magnetic Reynolds number Rm 
(that can be as high as Rm ^ 10^®) present a big problem for numerical simulations that cannot possibly get even 
close to the astrophysically- motivated numbers. The currently available 3D simulations can have Re and Rm up to 
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~ 10^. Both scale as the size of the box to the first power, while the computational effort increases as the fourth power 
(3 coordinates + time), so the brute force approach cannot begin to resolve the controversies related, for example, to 
ISM turbulence. 

We expect that observational studies of turbulence velocity spectra will provide important insights into ISM physics. 
Even in the case of much more simple oceanic (essentially incompressible) turbulence, studies of spectra allowed to 
identify meaningful energy injection scales. In interstellar, intra-cluster medium, in addition to that, we expect to see 
variations of the spectral index arising from the variations of the degree of compressibility, magnetization, interaction 
of different interstellar phases etc. 

How to get the turbulence spectra from observations is a problem of a long standing. While density fluctuations are 
readily available through both interstellar scincillations and studies of column density maps, the more coveted velocity 
spectra have been diSicult to obtain reliably until very recently. 

Turbulence is associated with fluctuating velocities that cause fluctuations in the Dopplcr shifts of emission and 
absorption lines. Observations provide integrals of either emissivities or opacities, both proportional to the local 
densities, at each velocity along the line of sight. It is far from trivial to determine the properties of the underlying 
turbulence from the observed spectral line. 

Centroids of velocity (Munch 1958) have been an accepted way of studying turbulence, although it was not clear 
to when and to what extend the measure really represents the velocity. Recent studies (Lazarian & Esquivel 2003, 
henceforth LE03, Esquivel & Lazarian 2005, Ossenkopf et al 2006, Esquivel ct al. 2007) have showed that the centroids 
are not a good measure for supersonic turbulence, which means that while the results obtained for HII regions (O'Dell 
& Castaneda 1987) are probably OK, those for molecular clouds are unreliable. 

An important progress in analytical description of the relation between the spectra of turbulent velocities and the 
observable spectra of fluctuations of spectral intensity was obtained in Lazarian & Pogosyan (2000, henceforth LPOO). 
This description paved way to two new techniques, which were later termed Velocity Channel Analysis (VGA) and 
Velocity Coordinate Spectrum (VCS). 

The techniques provide different ways of treating observational data in Position- Position- Velocity (PPV) data cubes. 
While VGA is based on the analysis of channel maps, which are the velocity slices of PPV cubes, the VCS analyses 
fluctuations along the velocity direction. If the slices have been used earlier for turbulence studies, although the 
relation between the spectrum of intensity fluctuations in the channel maps and the underlying turbulence spectrum 
was unknown, the analysis of the fluctuations along the velocity coordinate was initiated by the advent of the VCS 
theory. 

With the VCA and the VCS one can relate both observations and simulations to turbulence theory. For instance, 

the aforementioned turbulence indexes are very informative, e.g. velocity indexes steeper than the Kolmogorov value 
of —5/3 are likely to reflect formation of shocks, while shallower indexes may reflect scale-dependent suppression of 
cascading (see Beresnyak & Lazarian 2006 and ref. therein). By associating the variations of the index with different 
regions of ISM, e.g. with high or low star formation, one can get an important insight in the fundamental properties 
of ISM turbulence, its origin, evolution and dissipation. 

The absorption of the emitted radiation was a concern of the observational studies of turbulence from the very start 
of the work in the field (see discussion in Munch 1999). A quantitative study of the effects of the absorption was 
performed for the VCA in Lazarian & Pogosyan (2004, henceforth LP04) and for the VCS in Lazarian & Pogosyan 
(2006, henceforth LP06). In LP06 it was stressed that absorption lines themselves can be used to study turbulence. 
Indeed, the VCS is a unique technique that does not require a spatial coverage to study fluctuations. Therefore 
individual point sources sampling turbulent absorbing medium can be used to get the underlying turbulent spectra. 

However, LP06 discusses only the linear regime of absorption, i.e. when the absorption lines arc not saturated. This 
substantially limits the applicability of the technique. For instance, for many optical and UV absorption lines, e.g. Mg 
II, SII, Sill the measured spectra show saturation. This means that a part of the wealth of the unique data obtained 
e.g. by HST and other instruments cannot be handled with the LP06 technique. 

The goal of this paper is to improve this situation. In particular, in what follows, we develop a theoretical description 
that allows to relate the fluctuations of the absorption line profiles and the underlying velocity spectra in the saturated 
regime. 

Below, in §2 we describe the setting of the problem we address, while our main derivations are in §3. The discussion 
of the new technique of turbulence study is provided in §4, while the summary is in §5. 

2. ABSORPTION LINES AND MATHEMATICAL SETTING 

While in our earlier publications (LPOO, LP04, LP06) concentrated on emission line, in particular radio emission 
lines, e.g. HI and CO, absorption lines present the researchers with well defined advantages. For instance, they allow 
to test turbulence with a pencil beam, suffer less from uncertainties in path length. In fact, studies of absorption 

features in the spectra of stars have proven useful in outlining the gross features of gas kinematics in Milky way. 
Recent advances in sensitivity and spectral resolution of spectrographs allow studies of turbulent motions. 

Among the available techniques, VCS is the leading candidate to be used with absorption lines. Indeed, it is only 
with extended sources that the either centroid or VCA studies are possible. At the same time, VCS makes use not of 
the spatial, but frequency resolution. Thus, potentially, turbulence studies are possible if absorption along a single line 
is available. In reality, information along a few lines of sight, as it shown in Fig 1 is required to improve the statistical 
accuracy of the measured spectrum. Using the simulated data sets Chepurnov & Lazarian (2006ab) experimentally 
established that the acceptable number of lines ranges from 5 to 10. 
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For weak absorption, the absorption and emission lines can be analyzed in the same way, namely, the way suggested 
in LP06. For this case, the statistics to analyse is the squared Fourier transform of the Doppler-shifted spectral line, 
irrespectively of the fact whether this is an emission or an absorption spectral line. Such a "spectrum of spectrum" 
is not applicable for saturated spectral lines, which width is still determined by the Doppler broadening. It is known 
(see Spitzer 1978) that this regime corresponds to the optical depth r ranging from 10 to 10^. The present paper will 
concentrate on this regime^. 

Observatiqns in absorption line 




Cluster of 569ra 



TurftMiierii .^beeneiiOii ^fm uire in pure higin^soiuiiOTi iricae 

Fig. 1. — Absorption line study of turbulence in a cloud with VCS. 



Consider the problem in a more formal way. Intensity of the absorption line at frequency is given as 

/(z/, i^o) = /oe-"'"^''"' (1) 
where T(t/) is the optical depth. In the limit of vanishing intrinsic width of the line — i^o) = cw5(i^ — i^q), the 
frequency spread of r(z^, j/q) is determined solely by the Doppler shift of the absorption frequency from moving atoms. 
The number density of atoms along the line of sight moving at required velocity v ~ -^(y — vq) is 

= / dzp{z)4>{v - u{z) ~ Vreg) (2) 

Jo 

where (f) is the thermal distribution centered at every point z at the local mean velocity that is determined by the sum 
of turbulent and regular flow at that point. This is the density in PPV coordinate that we introduced in LPOO, so 
t{v) — a{i'o)ps{v). The intrinsic line width is accounted for by the convolution 

r(w) = a(i^o) y dwps{w)4>i{w - v) (3) 

or, in more detail, 

t{v) = a{va) j dzn{z) j dw(t){w - u{z) - Vgai{z))4>i{w - v) . (4) 
With intrinsic profile given by the Lorenz form 0i(u> — v)= („,_°,')2_|_„:i , the inner integral gives the shifted Voigt profile 
H{v - uix) - Vgai{x)) = / dw- -J— ^ exp — ^ , 5 

^ It is known that for r larger than 10^ the line width is determined by atomic constants and therefore it does not carry information 
about turbulence. 
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so we have another representation 

t{v) = a(fo) / dzn{x)H{v — u{x)) . (6) 

We clearly see from Eq. ([6]) that the line is affected both by Doppler shifts^ and atomic constants. 

3. FLUCTUATIONS OF OPTICAL DEPTH 
3.1. Statistics of Optical Depth Fluctuations 

The optical depth as a function of frequency contains fluctuating component arising from turbulent motions and 
associated density inhomogeneities of the absorbers. Statistics of optical depth fluctuations along the line of sight 
therefore carries information about turbulence in ISM. 

The optical depth is determined by the density of the absorbers in the PPV space, ps- In our previous work we have 
studied statistical properties of ps in the context of emission lines, using both structure function and power spectrum 
formalisms. Absorption lines demonstrate several important differences that warrant separate study. 

Firstly, our ability to recover the optical depth from the observed intensity 

/(i.) = /oe-"M + AT (7) 

depends on the magnitude of the absorption as well as sensitivity of the instrument and the level of measurement noise 
N. 

For lines with low optical depth tq < — ln(A^//o) we can in principle measure the optical depth throughout the whole 
line. At higher optical depths, the central part of the line is saturated below the noise level and the useful information 
is restricted to the wings of the line. This is the new regime that is the subject of this paper. 

In this regime the data is available over a window of frequencies limited to velocities high enough so that t{v) < 

— ln(7V//o) but not as high as to have Lorentz tail deflne the line. Higher the overall optical depth, narrower are the 
wings (following Spitzer, at tq > 10^ the wings are totally dominated by Lorentz factor). We shall denote this window 
by W{v — wqj ^) where uq is the velocity that the window is centered upon (describing frequency position of the wing) 
and A is the wing width. It acts as a mask on the "underlying" data 

t{v) ->t(w)H/(u-wo,A) . (8) 

Secondly, fluctuations in the wings of a line are superimposed on the frequency dependent wing profile. In other 
words, the statistical properties of the optical depth are inhomogeneous in this frequency range, with frequency 
dependent statistical mean value. While fluctuations of the optical depth 5t{v) that have origin in the turbulence can 
still be assumed to be statistically homogeneous, the mean profile of a wing must be accounted for. 

What statistical descriptors one should chose in case of line of sight velocity data given over limited window? Primary 
descriptors of a random field, here t{v)^ are the ensemble average product of the values of the field at separated points 

— the two point correlation function 

ir{vi,V2) ^ {t{vi)t{v2)) (9) 

and, reciprocally, the average square of the amplitudes of its (Fourier) harmonics decomposition — the power spectrum 

Pr{K,lK,2) = (r(fc„,i),T*(fc„,2)) (10) 

In practice these quantities are measurable if one can replace ensemble average by averaging over different positions 
which relies on some homogeneity properties of stochastic process. We assume that underlying turbulence is homoge- 
neous and isotropic. This does not make the optical depth to be statistically homogeneous in the wings of the line, 
but allows to introduce the fluctuations of r on the background of the mean profile t(w), At(w) — t{v) — f{v), which 
are (LP04) 

eAr(f) = (AT(z;i)Ar(w2)) , V = Vi - V2 (11) 
-PAr(fc„,lfc«,2) = (AT(fc,)AT*(fc„)) 5{K, - K,) (12) 

Homogeneous correlation function depends only on a point separation and amplitudes of distinct Fourier harmonics 
are independent. The obvious relations are 

(.t{vi,V2)^(,\t{v) +f{vi)f{v2) (13) 
Pr{kv^K,2)^PAAK)+T{Kl)f*{k^2) (14) 

Although mathematically the power spectrum is just a Fourier transform of the correlation function 

PAr(fc.) = j e'''^''Ur{v)dv , (15) 
which of them is best estimated from data depends on the properties of the signal and the data. 



^ Speaking formally, one can always make use of the known atomic constants and get inside into turbulence. This is not practical, 
however, for an actual spectrum in the presence of noise. 
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The power spectrum carries information which is locahzed to a particular scale and as such is insensitive to processes 
that contribute outside the range of scales of interest, in particular to long-range smooth variations. On the other hand, 
determination of Fourier harmonics is non-local in configuration space and is sensitive to specifics of data sampling - 
the finite window, discretization, that all lead to aliasing of power from one scales to another. The issue is severe if 
the aliased power is large. 

Conversely, the correlation function is localized in configuration space and can be measured for non-uniformly 
sampled data. However, at each separation it contains contribution from all scales and may mix together the physical 
effects from different scales. In particular, ^{v) is not even defined for power law spectra P{ky) ~ with index 
n > 1 (for one dimensional data) This limitation is relieved if one uses the structure function 

DaAv) ^ {{At{vi) - At{v2))^) (16) 

instead, which is well defined for n < 3. The structure function can be thought of as regularized version of the 
correlation function 

DAr{v)^2{UriO)-Ur{v)) (17) 

that is related to the power spectrum in the same way as the correlation function, if one excludes the fc^, = mode. 
Velocity Coordinate Spectrum studies of LP06 demonstrated that the expected one dimensional spectrum of PPV 

density fluctuations along velocity coordinate that arise from turbulent motions is P{kv) ~ where m is the index 

of line-of-sight component of the velocity structure function. For Kolmogorov turbulence m = 2/3 and for turbulent 
motions dominated by shocks m = 1/2. These spectra are steep ~ k~'^ , ~ fc""* which makes the direct measurement of 
the structure functions impractical (although for m > 2/3 the structure function can be defined). At the same time, 
in our present studies we deal with a limited range of data in the wings of the absorption lines, which complicates the 
direct measurements of the power spectrum. 

Below we first describe the properties of the power spectrum PArik-u) in this case, and next develop the formalism 
of higher order structure functions. 

3.2. Power Spectrum of Optical Depth Fluctuations 

Let us derive the power spectrum of the optical depth fluctuations, PT{ky) = (T(fc„)T*(fc„)) Here ky is wave number 
reciprocal to the velocity (frequency) separation between two points on the line-of-sight and angular brackets denote 
an ensemble averaging. 

Fourier transform of the eq. ^ with respect to velocity is 

r(fc„) = a(i/o) / dzn{x) [ dkle-^^>e-^e-<''^^'>W{ky-k'^,vo,A) (18) 
and the power spectrum 







Prikv) = (^a{i^o)^ J dzi J dz2n{zi)n{z2)x (19) 
dkl / dfc:^e-<l^^l+l''"l^°e- 2" ^' e-'iKu(zi)-Kuiz,)]-^y^f^^ _ k'jw*ik, - k'^) \ 

which is useful to express using average velocity u+ = {u{zi) + u{z2))/2 and velocity difference u = u{zi) — u{z2), as 
well as correspondent variables for the wave numbers fc+ = (fc(, + fc,")/2 and k~ = k'^ — fc", as 



P^{k^)^[a{vaf I dzi dz2n{zi)n{z2) I dfc+e-'=^ '^e''^^" (20) 
"'0 



The fluctuating, random quantities, over which the averaging is performed are the density n{z) and the line-of-sight 
component of the velocity of the absorbers u{z), varying along the line of sight . In our earlier papers (see LPOO, 
LP04) we argued that in many important cases they can be considered as uncorrelated between themselves, so that 

(n(zi)n(^2)e''=-^"e^'=^' "+) = ^{\z, - z2|)e-^'=:'^^-(l-i--2l)g-ife-^[i5.(S)-D.(|.i-..|)/2] ^ (31) 

where ^(|zi — Z2I) is the correlation function of the density of the absorbers and Dz{\zi — 2^21) is the structure function 
of their line-of-sight velocity due to turbulent motions. Dz is expected to saturate at the value Ds{S) for separations 

^ In physically realistic situations the power law will not, of course, extend infinitely to large scales. Mathematical divergence of the 
correlation function in practice mean that for steep spectra the largest present scales the correlation at all, even small separations v. 

^ note, that due to having measurements in the finite window, the amplitudes at different waves numbers are, in general, correlated, 
(r(fc„)r*(fc(,)) 7^ 0, ky ^ Here we restrict ourselves to diagonal terms only. 
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of the size of the absorbing cloud S. The dependence of ^(z) and Dz{z) only on spatial separation between a pair of 
absorbers reflects the assumed statistical homogeneity of the turbulence model. Introducing z — zi~Z2 and performing 
integration over z+ = (zi + Z2)/2 one obtains ^ 



s 

dz{S 



dk, 



+ correction terms {see Appendix) 

2/3, = Dz{S) — Dz{\zi — Z2|)/2 + /9 and symmetrized window W is defined in the 



(22) 



where = Dz{\zi — Z2I 
Appendix. 

If one has the whole line available for analysis, the masking window will be flat with (5(fc)-function like Fourier 
transform. The combination of the windows in the power spectrum will translate to d{ky ~ fc+)(5(fc~) and 



dz[S-\z\)i{z)e 



(23) 



Masking the data has the effect of aliasing modes of the large scales that exceed the available data range, to shorter 
wavelength. This is represented by the convolution with Fourier image of the mask. Secondary effect is the contribution 
of the modes with different wave numbers /c^ 7^ to the diagonal part of the power spectrum. This again reflects the 
situation that different Fourier components are correlated in the presence of the mask. 

To illustrate the effects of the mask, let us assume that we select the line wing with the help of a Gaussian mask 
centered in the middle of the wing at v\ 



W(v- 



2A^ 



-^;i,A): 

Ty(fc„) = %/2^Ae-*'=^'"i( 



that gives 



All integrals can then be carried out to obtain 



W = 47rA2 cos -^vx 



X exp 



dz{S-z)£,{z) 



' exp 



,-(fc„-fe+)2A2 



2a^ 



D- 



2A2 



exp 



+ 2A2)(D- 



D- 



2A2 



Erfc 



2A2) 
^/2{a 



exp 



-(fe„+/£+)^A = 



The following limits (taking k^ > 0) are notable: 
A ^ 00 : 



P{k,)oca{i^ofS 
a -> : 



dz{S — z)^{z) exp 



y/D- + 2A2 



2A2 



(24) 
(25) 

(26) 



(27) 



(28) 



(29) 



P{ky)<xa{iyofS dz{S z)^{z) 



A^ exp 



^(L>+ + 2A2)(L>- +2A2 



: exp 



D- 



2A2 



The last expression particularly clearly demonstrates the effect of the window, which width in case of the line wing 
is necessarily necessarily limited by A2 < D{S) + 2(3. The power spectrum is corrupted at scales fc„ < A^^, but 
still maintains information about turbulence statistics for fc„ 3> A~^. Indeed, in our integral representation the power 

spectrum at k^ is determined by the linear scales such that -j^^^EJ^v ^ which translates into D~ < 2A'^/ {k'^A^ — 1). 
Thus, if fc„A 3> 1 over all scales deflning power at k^ one has A ^ and there is no significant power aliasing.^ 
For intermediate scales there is a power aliasing as numerical results demonstrate in Figure [2] Doppler broadening 
described by D~ incorporates both turbulent and thermal effects. Thermal effects are especially important in case of 
narrow line wings, since the range of the wavenumbers relatively unaffected by both thermal motions and the mask 
is limited l/sqrt/3 > ky > A^^ and exists only for relatively wide wings A > ^/p. For narrower wings the combined 
turbulent and thermal profile must be fitted to the data, possibly determining the temperature of the absorbers at 

^ Absolute values in the Fourier image of the Lorentz transform require separate consideration of each quadrant of wave numbers. This 
is carried out in the Appendix. 

® In this regime the first factor is essentially constant since Z)+ (2:) = (S) + 2(3 — varies little, (2:) Z)+(0) in the interval 

2/3 < < 2fc^^ < 2A2 < D(S) + 2(3. 
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1 2 5 10 20 50 100 1 2 5 10 20 50 100 



Fig. 2. — Power spectrum of the optical depth fluctuations from line wings. All parameters, ky, A, ''I a^r^ dimensionless, in the 
units of y/D~z(S), the variance of the turbulent velocity at the scale of the cloud. Intrinsic line broadening is neglected. Only the effect 
of the turbulent motions and not spatial inhomogeneity of the absorbers is taken into account. The underlying scaling of the turbulent 
velocities is Kolmogorov, m = 2/3. The left panel illustrates the power aliasing due to finite width of the window. The power spectrum 
is plotted, from top to bottom, for A = 1, 0.2.0.1, i.e the widths of the wing ranges from the complete line to one-tenth of the line width. 
The straight line shows the power law P{ky) oc fc„ '^^"^ = k^^ expected under ideal observational circumstances. One finds that for the 
ideal Gaussian mask the underlying spectrum is recovered for k^ > 3A~^. The right panel shows the modification of the spectrum due to 
thermal broadening, which is taken at the level = ^ \/ Dz{S). Thermal effects must be accounted for for k^ > l/(3\//3)- 

the same time. This recipe is hmited by the assumption that the temperature of the gas is relatively constant for the 
absorbers of a given type. 

We should note that the Gaussian window provides one of the ideal cases, limiting the extend of power aliasing since 
the window Fourier image falls off quickly. One of the worst scenarios is represented by sharp top hat mask, which 
Fourier image falls off only as k~'^ spreading the power from large scales further into short scales. For steep spectra 
that we have in VCS studies all scales may experience some aliasing. This argues for extra care while treating the line 
wings through power spectrum or for use of alternative approaches. 



3.3. Second order structure function 

Second order structure function provides an alternative to power spectrum measurement in case of steep spectra 
with the data limited to the section of the lines. The second order structure function of the fluctuations of the optical 
depth can be defined as 



1 



dDAr{v) = - ( (At(wi +v) + At{vi 



v) - 2At{vi)Y 

It represents additional regularization of the correlation function beyond the ordinary structure function 

dDAr{v) = S^Ar (0) - i^Av) + ^Ar{2v) 



(30) 

(31) 
(32) 



Dat{v) is proportional to three dimensional velocity space density structure function at zero angular separations 
dp(0, v), discussed in LP06. Using the results of LP06 for dp(0, v) we obtain 



dDAr{v)-r-^ 



S 



dz 



1 



£|7+m/2 



3 — 4 cxp 



1 



D^{S)m \S 



2\z\ 



- cxp 

24p— 6 



p 



2v' 



(33) 



where p = (1 — j)/m — 1/2 > and 7 is the correlation index that describes spatial inhomogeneities of the absorbers. 

1/2 

To shorten intermediate formulas, the dimensionless quantities v = v/ Dz {S), z = z/S, f = r/S are introduced. The 
first term in the expansion contains information about the underlying field, while the power law series represent the 
effect of boundary conditions at the cloud scale. In contrast to ordinary correlation function they are not dominant 
until p > 2, i.e for (1 — j)/m < 5/2 the second order structure function is well defined. When turbulent motions 
provide the dominant contribution to optical depth fluctuations, 7 = 0, we see that measuring the dD^. one can recover 
the turbulence scaling index if m > 2/5, which includes both interesting cases of Kolmogorov turbulence and shock- 
dominated motions. This condition is replaced by m > 2/5(1 — 7) if the density fluctuations, described by correlation 
index 7, are dominant. 

At sufficiently small scales the second order structure function has the same scaling as the first order one 



dDAriv) oc 



Dz{S)m \S 



2P- 



2^-p) r(-p) . v^p 



(34) 
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A practical issue of measuring the structure functions directly in the wing of the line is to take into account the line 
profile. The directly accessible 

dDrivuv) = ^ ((r(vi +v)+ t{vi -v)- 2T{vi)f) (35) 
is related to the structure function of the fluctuations as 

dDAr(v) = dDr{vi,v) - i {f{vi +v)+ f{vi -v)- 2f{vi)f (36) 

where the mean profile of the optical depth f is related to the mean profile of PPV density ps given in the Appendix 
B of LP06. At small separations v, the correction to the structure function due to mean profile behaves as and is 
subdominant. 

The price one pays when utilizing higher-level structure function is their higher sensitivity to the noise in the data. 
While correlation function itself is not biased by the noise except at at zero separations (assuming noise is uncorrelatcd 

([t(z;i) + N{vi)][T{v, +v) + N{v, + v)]) = ^v) + {N^)6{v) (37) 
already the structure function is biased by the noise which contributes to all separations 

([t(wi) + N{vi) - t(wi +v)- N{vi + w)]2) = Dr{v) + 2{N^) (38) 
This effect is further amplified for the second order structure function 

^ ^[r(«i +v) + t{vi - u) - 2t{vi) + N{vi + v) + N{vi - v) - 2N{vi)]^'^ = dDr{v) + 3{N^) (39) 

The error in the determination of the structures function of higher order due to noise also increases. 

3.4. Comparison of the approaches 

Structure functions and power spectra are used interchangeably in the theory of turbulence (see Monin & Yaglom 
1975). However, complications arise when spectra are "extremely steep", i.e. the corresponding structure function of 
fluctuation grows as .x^, H > 2. For such random fields, one cannot use ordinary structure functions, while the one 
dimensional Fourier transforms that is employed in VCS corresponds to the power spectrum of k~^^^ is well defined. 

As a rule, one docs not have to deal with so steep spectra in theory of turbulence (see, however, Cho et al. 2002 and 
Cho & Lazarian 2004). Within the VCS, such "extremely steep" spectra emerge naturally, even when the turbulence 
is close to being Kolmogorov. This was noted in LP06, where the spectral approach was presented as the correct one 
to studying turbulence using fluctuations of intensity along v-coordinatc. 

The disadvantage of the spectral approach is when the data is being limited by a non-Gaussian window function. 
Then the contributions from the scales determined by the window function may interfere in the obtained spectrum at 
large ky. An introduction of an additional more narrow Gaussian window function may mitigate the effect, but limits 
the range of ky for which turbulence can be studied. Thus, higher order structure functions (see the subsection above), 
is advantageous for the practical data handling. 

In terms of the VGA theory, we used mostly spectral description in LPOO, while in LPG4, dealing with absorption, 
we found advantageous to deal with real rather than Fourier space. In doing so, however, we faced the steepness of 
the spectrum along the v-coordinate and provided a transition to the Fourier description to avoid the problems with 
the "extremely steep" spectrum. Naturally, our approach of higher order structure functions is applicable to dealing 
with the absorption within the VGA technique. 

4. DISCUSSION 
4.1. Simplifying assumptions 

In the paper above we have discussed the application of VGS to strong absorption lines. The following assumption 
were used. First of all, considering the radiative transfer we neglected the effects of stimulated emission. This 
assumption is well satisfied for optical or UV absorption lines (see Spitzer 1979). Then, we assumed that the radiation 
is coming from a point source, which is an excellent approximation for the absorption of light of a star or a quasar. 
Moreover, we disregarded the variations of temperature in the medium. 

Within our approach the last assumption may be most questionable. Indeed, it is known that the variations of 
temperature do affect absorption lines. Nevertheless, our present study, as well as our earlier studies, prove that the 
effects of the variations of density are limited. It is easy to see that the temperature variations can be combined 
together with the density ones to get effective renormalized "density" which effects we have already quantified. 

Our formalism can also be generalized to include a more sophisticated radiative transfer and the spatial extend of 
the radiation source. In the latter case we shall have to consider both the case of a narrow and a broad telescope 
beam, the way it has been done in LP06. Naturally, the expressions in LPG6 for a broad beam observations can be 
straightforwardly applied to the absorption lines, substituting the optical depth variations instead of intensities. The 
advantage of the extended source is that not only VGS, but also VGA can be used (see Deshpande et al. 2000). As a 
disadvantage of an extended source is the steepening of the observed v-coordinate spectrum for studies of unresolved 
turbulence. This, for instance, may require employing even higher order structure functions, if one has to deal with 
windows arising from saturation of the absorption line. 
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In LP06 wc have studied the VCS technique in the presence of absorption and formulated the criterion for the 
fluctuations of intensity to reliably reflect the fluctuations in turbulent velocities. In this paper, however, we used the 
logarithms of intensities and showed that this allows turbulence studies beyond the regime, at which fluctuations of 
intensity would be useful. 

A similar approach, potentially, is applicable to the VCS for studies of emitted radiation in the presence of significant 
absorption. Indeed, the equation of the radiative transfer provides us in this case with 

Jx(«) = - [l -c-"''^^^'^)] , (40) 
a L J 

where the first term is a constant, that can be subtracted. If this is done, the logarithm can be taken of the intensity. 
Potentially, this should allow VCS studies of emission lines where, otherwise, absorption distorts the statistics. 

The difficulty of such an approach is the uncertainty of the base level of the signal. Taking logarithm is a non-linear 
operation that may distort the result, if the base level of signal is not accounted for properly. However, the advantage 
of the approach that potentially it allows studies of velocity turbulence, when the traditional VGA and VCS fail. 
Further research should clarify the utility of this approach. 

4.2. Prospects of Studying ISM Turbulence with Absorption Lines 

The study of turbulence using the modified VCS technique above should be reliable for optical depth r up to f 0'^. 
For this range of optical depth, the line width is determined by Doppler shifts rather than the atomic constants. While 
formally the entire line profile provides information about the turbulence, in reality, the flat saturated part of the 
profile will contain only noise and will not be useful for any statistical study. Thus, the wings of the lines will contain 
signal. 

As several absorption lines can be available along the same line of sight, this allows to extend the reliability of 
measurements combining them together. We believe that piecewise analyses of the wings belonging to different 
absorption lines is advantageous. 

The actual data analysis may employ fitting the data with models, that, apart from the spectral index, specify the 
turbulence injection scale and velocity dispersion, as this is done in Chepurnov et al. (2006). 

Note, that measurements of turbulence in the same volume using different absorption lines can provide comple- 
mentary information. Formally, if lines with weak absorption, i.e. tq < I are available, there is no need for other 
measurements. However, in the presence of inevitable noise, the situation may be far from trivial. Naturally, noise of a 
constant level, e.g. instrumental noise, will affect more weak absorption lines. The strong absorption lines, in terms of 
VCS sample turbulence only for sufficiently large This limits the range of turbulent scales that can be sampled with 
the technique. However, the contrast that is obtained with the strong absorption lines is higher, which provides an 
opportunity of increasing signal to noise ratio for the range of fc„ that is sampled by the absorption lines. If, however, 
a single strong absorption line is used, an analogy with a two dish radio interferometer is appropriate. Every dish of 
the radio interferometer samples spatial frequencies in the range approximately [1/A, 1/d], where A is the operational 
wavelength, d is the diameter of the dish. In addition, the radio interferometer samples the spatial frequency 1/D, 
where D is the distance between the dishes. Similarly, a strong absorption line provides with the information on 
turbulent velocity at the largest spatial scale of the emitting objects, as well as the fluctuation corresponding to the 

scales knjindow ; ^abs • 

In LP06 we concentrated on obtaining asymptotic regimes for studying turbulence. At the same time in Chepurnov 
et al. (2006) fitting models of turbulence to the data was attempted. In the latter approach non-power law observed 
spectra can be used, which is advantageous for actual data, for which the range of scales in ky is rather limited. Indeed, 
for HI with the injection velocities of 10 km/s and the thermal velocities of I km/s provides an order of magnitude 
of effective "inertial range". Correcting for thermal velocities one can increase this range by a factor, which depends 
on the signal to noise ratio of the data. Using heavier species rather than hydrogen one can increase the range by a 
factor {m heavy /i^Hy^"^- This may or may not be enough for observing good asymptotics. 

We have seen in ?? that for absorption lines the introduction of windows determined by the width of the line wings 
introduces additional distortions of the power spectrum. However, this is not a problem if, instead of asymptotics, 
fitting of the model is used. Compared to the models used in Chepurnov et al. (2006) the models for absorption lines 
should also have to model the window induced by the absorption. The advantage is, however, that absorption lines 
provide a pure pencil beam observations. 

4.3. Comparison and synergy with other techniques 

Formally, there exists an extensive list of different tools to study turbulence that predated our studies (see Lazarian 
1999 and ref. therein). However, a closer examination shows that this list is not as impressive as it looks. Moreover, our 

research showed that some techniques may provide confusing, if not erroneous, output, unless theoretical understanding 
of what they mcasiire is achieved. For instance, we mentioned in the introduction an example of the erroneous 
application of velocity ccntroids to supersonic molecular cloiid data. Note, that clumps and shell finding algorithms 
would find a hierarchy of clumps/shells for synthetic observations obtained with incompressible simulations. This calls 
for a more cautious approach to the interpretation of the results of some of the accepted techniques. 

For instance, the use of different wavelets for the analysis of data is frequently treated in the literature as different 
statistical techniques of turbulence studies (Gill & Henriksen 1990, Stutzki et al. 1998, Cambresy 1999, Khahl et al. 
2006), which creates an illusion of an excessive wealth of tools and approaches. In reality, while Fourier transforms 
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use harmonics of e'^'' , wavelets use more sophisticated basis functions, which may be more appropriate for problems 
at hand. In our studies we also use wavelets both to analyze the results of computations (see Kowal & Lazarian 
2006a) and synthetic maps (Ossenkopf et al. 2006, Esquivel et al. 2007), along with or instead of Fourier transforms 
or correlation functions. Wavelets may reduce the noise arising from inhomogeneity of data, but we found in the 
situations when correlation functions of centroids that we studied were failing as the Mach number was increasing, a 
popular wavelet (A-variance) was also failing (cp. Esquivel & Lazarian 2005, Ossenkopf et al. 2006, Esquivel et al. 
2007). 

While in wavelets the basis functions are fixed, a more sophisticated technique. Principal Component Analysis 
(PCA), chooses basis functions that are, in some sense, the most descriptive. Nevertheless, the empirical relations 
obtained with PCA for extracting velocity statistics provide, according to Padoan et al. (2006), an uncertainty of the 
velocity spectral index of the order 0.5 (see also Brunt et al. 2003), which is too large for testing most of the turbulence 
theories. In addition, while our research in LPOO shows that for density spectra Ep ~ for a < 1 both velocity 
and density fluctuations influence the statistics of PPV cubes, no dependencies of PPV statistics on density have been 
reported so far in PCA studies. This also may reflect the problem of finding the underlying relations empirically with 
data cubes of limited resolution. The latter provides a special kind of shot noise, which is discussed in a number of 
papers (Lazarian et al. 2001, Esquivel et al. 2003, Chepurnov & Lazarian 2006a). 

Spectral Correlation Function (SCF) (sec Rosolowsky ct al. 1999 for its original form) is another way to study 
turbulence. Further development of the SCF technique in Padoan ct al. (2001) removed the adjustable parameters 
from the original expression for the SCF and made the technique rather similar to VCA in terms of the observational 
data analysis. Indeed, both SCF and VCA measure correlations of intensity in PPV "slices" (channel maps with a given 
velocity window Au), but if SCF treats the outcome empirically, the analytical relations in Lazarian & Pogosyan (2000) 
relate the VCA measures to the underlying velocity and density statistics. Mathematically, SCF contains additional 
square roots and normalizations compared to the VCA expressions. Those make the analytical treatment, which 
is possible for simpler VCA expressions, prohibitive. One might speculate that, similar to the case of conventional 
centroids and not normalized centroids introduced in Lazarian & Esquivel (2003), the actual difference between the 
statistics measured by the VCA and SCF is not significant. 

In fact, we predicted several physically-motivated regimes for VCA studies. For instance, slices arc "thick" for 
eddies with velocity ranges less than Au and "thin" otherwise. VCA relates the spectral index of intensity fluctuations 
within channel maps to the thickness of the velocity channel and to the underlying velocity and density in the emitting 
turbulent volume^. In the VCA these variations of indexes with the thickness of PPV "slice" are used to disentangle 
velocity and density contributions. We suspect that similar "thick" and "thin" slice regimes should be present in the 
SCF analysis of data, but they have not been reported yet. While the VCA can be used for all the purposes the SCF 
is used (e.g. for an empirical comparisons of simulations and observations), the opposite is not true. In fact, Padoan 
et al. (2004) stressed that VCA eliminates errors inevitable for empirical attempts to calibrate PPV fluctuations in 
terms of the underlying 3D velocity spectrum. 

VCS is a statistical tool that uses the information of fluctuations along the velocity axis of the PPV. Among all 
the tools that use spectral data, including the VCA, it is unique, as it does not require spatial resolution. This is 
why, dealing with the absorption lines, where good spatial coverage is problematic, we employed the VCS. Potentially, 
having many sources sampling the object one can create PPV cubes and also apply the VCA technique. However, this 
requires very extended data sets, while for the VCS sampling with 5 or 10 sources can be sufficient for obtaining good 
statistics (Chepurnov & Lazarian 2006a). 

We feel that dealing with the ISM turbulence, it is synergctic to combine different approaches. For the wavelets used 
their relation with the underlying Fourier spectrum is usually well defined. Therefore the formulation of the theory 
(presented in this work, as well as, in our earlier papers in terms of the Fourier transforms) in terms of wavelets is 
straightforward. At the same time, the analysis of data with the wavelets may be advantageous, especially, in the 
situations when one has to deal with window functions. 

5. SUMMARY 

In the paper above we have shown that 

• Studies of turbulence with absorption lines are possible with the VCS technique if, instead of intensity /(z^), one 
uses the logarithm of the absorbed intensity log/Q5s(i^), which is equivalent to the optical depth t{v). 

• In the weak absorption regime, i.e. when the optical depth at the middle of the absorption line is less than unity, 
the analysis of the t{v) coincides with the analysis of intensities of emission for ideal resolution that we discussed 
in LP06. 

• In the intermediate absorption retime, i.e. when the optical depth at the middle of the absorption line is larger 
than unity, but less than 10^, the wings of the absorption line can be used for the analysis. The saturated part 
of the line is expected to be noise dominated. 

• The higher the absorption, the less the portion of the spectrum corresponds to the wings available for the 

We showed that much of the eaxUer confusion stemmed from different observational groups having used velocity channels of different 
thicknesses (compare, e.g. , Green 1993 and Stanimirovic et al. 1999). 
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analysis. In terms of the mathematical setting this introduces and additional window in the expressions for the 
VCS analysis. However, the contrast of the small scale fluctuations increases with the decrease of the window. 

• For strong absorption regime, the broadening is determined by Lorentzian wings of the line and therefore no 
information on turbulence is available. 



APPENDIX 

DERIVATION OF THE POWER SPECTRUM Pr 
Following eqns. (j20|2ip the power spectrum of the optical depth is 

-Pr(fci;) = a(t'o)^ / dzi / dz2 ^{\zi - Z2\) I dk[, I dk'^ 
Jo Jo 



— oo ^ — oo 
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where = D,{\zi - zai) + 2/3 and D+ = D,{S) - D,{\zi - Z2\)/2 + /3 while k+ = (fc(, + k'^)/2 and fc- = k'^ - k'^. 
Since the mask is real, W*{kv) = W{—ki,). 

To deal with absolute values in the Lorentz transform, we split integration regions in quadrants I - > 0, A:" > 0), 
II - (fc(, < 0, A:" > 0), HI ~ < 0, fc" < 0) and IV - > 0, A:" < 0). Integration over quadrants III and IV can be 



folded into integration over regions I and II respectively by substitution fc(, — fc" — )■ — fc", k^ — )■ —k^,k^ 
Writing out only integration over ky 
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Changing variables of integration to fc+ and k^ 
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At the end, the integrals can be combined into the main contribution and the correction that manifests itself only 
when Lorentz broadening is significant. 
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where symmetrized window is 
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The final expression for Pr{ky) is then 

P^{k^) = a^S dz{S - \z\)^{z) dk+e 
Jq Jo 

/•S pcc 2 

+ a^S / |z|)C(z) / dfc+e-5fe^ d 



dk„ e 
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(A9) 
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